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Abstract. We consider several models of the damped oscillators in nonrelativistic quantum me- 
chanics in a framework of a general approach to the dynamics of the time-dependent Schrodinger 
equation with variable quadratic Hamiltonians. The Green functions are explicitly found in terms of 
elementary functions and the corresponding gauge transformations are discussed. The factorization 
technique is applied to the case of a shifted harmonic oscillator. The time-evolution of the expec- 
tation values of the energy related operators is determined for two models of the quantum damped 
oscillators under consideration. The classical equations of motion for the damped oscillations are 
derived for the corresponding expectation values of the position operator. 



1. An Introduction 

We continue an investigation of the one-dimensional Schrodinger equations with variable qua- 
dratic Hamiltonians of the form 

= -« (^) ^ + ^ w - ' (^) + ^ (^) ^) ' (1-1) 

where a (t) , & (t) , c (t) , and d{t) are real- valued functions of time t only; see Refs. [8J, |9], [22], 
[23], [25], [M], [35], and [36j for a general approach and currently known explicit solutions. Here 
we discuss elementary cases related to the models of damped oscillators. The corresponding Green 
functions, or Feynman's propagators, can be found as follows [8], [35] : 

^ = G{x,y,t) = , ^. ^r{a(t).^+m-y+,{t)y-) ^ (1^2) 
A/27ri/i (t) 



where 



1 d{t) 

a{t)= - TT^, (1-3) 



4a (t) /i (t) 2a (t) ' 
h{t) 



/? (t) = h (t) = exp ( - / (c (r) - 2d (r)) dr ) , (1.4) 







fiWA'W 2o(0) i„ {ij.'{T)f ^ 
and the function /i (t) satisfies the characteristic equation 

fx" - T (t) fi' + Aa (t) fx = (1.6) 
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with 

r{t) = -~2c + Ad, a{t)=ab-cd + d'^ + ^(--^] (1.7) 
subject to the initial data 

/i(0) = 0, /i' (0) = 2a (0) ^ 0. (1.8) 

More details can be found in Refs. [8] and [35] . The corresponding Hamiltonian structure is discussed 
in Ref. 0. 

The simple harmonic oscillator is of interest in many advanced quantum problems [16], [21], [26], 
and [32]. The forced harmonic oscillator was originally considered by Richard Feynman in his path 
integrals approach to the nonrelativistic quantum mechanics [12], [13], [S], [IS], and [IS]; see also 
[23] • Its special and limiting cases were discussed by many authors; see Refs. [S], [IZ], [IH], [21], 
[26], [38] for the simple harmonic oscillator and Refs. [1], [7], [18], [2Z], [31] for the particle in a 
constant external field and references therein. 

The damped oscillations have been analyzed to a great extent in classical mechanics; see, for 
example, Refs. [S] and [20]. In the present paper we consider the time-dependent Schrodinger 
equation 

= Hi; (1.9) 

with the following nonself-adjoint Hamiltonians 

H =^[p^ + x^) - \px (1.10) 

and 

H='^{p' + x')-\xp, (1.11) 
where p = —id/ dx, as quantum analogs of the damped oscillator. A related self-adjoint Hamiltonian 

H ='^(p^ + x^) -^{px + xp) (1.12) 

is also analyzed. Although discussion of a quantum damped oscillator is usually missing in the 
standard classical textbooks [21], [2S], and [32] among others, we believe that the models presented 
here have a significant value from the pedagogical and mathematical points of view. For instance, 
one of these models was crucial for our understanding of a "hidden" symmetry of the quadratic 
propagators in Ref. [9]. Moreover, our models show that fundamentals of quantum mechanics, 
such as evolution of the expectation values of operators and Ehrenfest's theorem, can be extended 
to the case of nonself-adjoint Hamiltonians. This provides, in our opinion, a somewhat better 
understanding of the mathematical foundations of quantum mechanics and can be used in the 
classroom. 

The paper is organized as follows. In section 2 we derive the propagators for the models of the 
damped oscillator (11.101) and (11. lip following the method of Ref. [8]. The corresponding gauge 
transformations are discussed in section 3. The next section is concerned with the separation of 
the variables for related model of a "shifted" linear harmonic oscillator fll.l2p . The factorization 
technique is applied to this oscillator in section 5. The time evolution of the expectation values of 
the energy related operators is determined for these quantum damped oscillators in section 6. The 
classical equations for the damped oscillations are derived for the expectation values of the position 
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operator in the next section. One more model of the damped oscillator with a variable quadratic 
Hamiltonian is introduced in section 8. The last section contains some remarks on the momentum 
representation. 



2. The First Two Models 
For the time- dependent Schrodinger equation: 

with a = b = ujq/2 and c = d = —A, the characteristic equation takes the form of the classical 
equation of motion for the damped oscillator [5] , [20] : 

H" + 2Xfi' + ulfi = 0, (2.2) 

whose suitable solution is as follows 



■^0 -Xt ■ 

— e SI 

u 

The corresponding propagator is given by 



/i = — e"^* sin ut, u = \/ul - > 0. (2.3) 

UJ 



G{x,y,t) = \ exp [[x + y ) cosut — Zxy) 

y zmujQsmujt \ 2ct;o sm cut 

xexp(^{x'-y')). (2.4) 



Indeed, directly from (ra-l lLi]) : 

cos cijt + A sin cut u 

The integral in (11.51) can be evaluated with the help of a familiar antiderivative 

dt sin t 



I 



+ C. (2.6) 



{Acost + Bsintf A{Acost + B sint) 
It gives 

UJ cos cut — X sin ut 
7W = 5 (2.7) 

with the help of the following identity 

uj"^ — Lul sin^ Lut = uJ^ cos^ u)t — A^ sin^ uot (2-8) 

and the propagator (12.41) is verified. A "hidden" symmetry of this propagator is discussed in Ref. [9]. 

The time-evolution of the squared norm of the wave function is given by 

/■oo 

||^(a;,t)f = / |^(x,t)|^ dx = e^*||^(x,0)|r (2.9) 



It is derived in section 6 among other things. We have discussed here the case tUg > A^. Two more 
cases, when = A^ and cUq < A^, are similar and the details are left to the reader. 
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In a similar fashion, the time-dependent Schrodinger equation of the form 

with a = h = ujq/2 and c = —A, ci = 0, has the characteristic equation 

/x" - 2A/i' + cu^/i = (2.11) 

with the solution 

;x = — e^*sincut, uj = J ujI - > 0. (2.12) 
The corresponding propagator is given by 



^ / N one f // 2 2\ \ 

G{x,y,tj = \ exp ((a; + y ) cosut — 2xy) 

y zmuosmut \2LJosmut 

xe^p(^{x'-y')] (2.13) 



,2uJo 

and the evolution of the squared norm is as follows 

\\i;{x,t)f = e-''U{x,0)f. (2.14) 

The solution of the Cauchy initial value problem 

^^ = H^, ^{x,0) = x{x) (2.15) 
for our models (12. ip and (I2.10p is given by the superposition principle in an integral form 

/oo 
G(x,y,t) x{y) dy (2.16) 
-oo 

for a suitable initial function % on i?; a rigorous proof is given in Ref. [35] . 

3. The Gauge Transformations 

The time-dependent Schrodinger equation 

= [^{p- Af + U + (p- A)V + W (p- A)')iP, (3.1) 

where p = i~^d/dx is the linear momentum operator and A = A{x,t) ^ U = U {x,t) , V = V {x,t) , 
W = W {x, t) are real- valued functions, with the help of the gauge transformation 

^ = e-^^(^'*V' (3.2) 

can be transformed into a similar form 

= (y (P - + U' + {p- A) V + W'{p- A)) ^' (3.3) 

with the new vector and scalar potentials given by 

A' = A + ^, U' = U-%, V' = V, W' = W. (3.4) 
ox at 

Here we consider the one-dimensional case only and may think of / as being an arbitrary complex- 
valued differentiable function. Also, the Hamiltonian in the right hand side of equation (13. ip is 
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not assumed to be self-adjoint. See Refs. |21j and [26] for discussion of the traditional case, when 
V =W = Q. 

An interesting special case of the gauge transformation related to this paper is given by 

^ = 0, ^/ = Y^^ v = -\x, iy = o, / = ^, (3.5) 

A' = 0, [/' = i^x2-^, V' = -\x, l^' = 0, (3.6) 

when the new Hamiltonian is 

H' = ^{p-A'f + U' + pV (3.7) 

+ 2x— + 1 



2 V <9a;2 / 2 \ dx 

and equation (12. ip takes the form 



The corresponding Green function is given by 

G(x,y,t) = A — exp ( — ((x'^ + y'^)cosujt — 2xy) 

Y 2TnuJosmujt \2uJosmujt 

xexp(^|^(x2-y2)^, u;=^Jcul-\^>0 (3.9) 

and the norm of the wave function is conserved with time. This can be established once again 
directly from our equations fll.2l) - fll.8l) . We leave the details to the reader. A traditional method 
of separation of the variables and using the Mehler formula for Hermite polynomials is discussed in 
the next section. The factorization technique is applied to this Hamiltonian in section 5. 

Equation (13.81) . in turn, admits another local gauge transformation: 

A = 0, f/=^x^ V = W = -^, / = -^, (3.10) 

A' = --, U' = ^x\ V' = W' = -^ (3.11) 
ujq 2 2 

and the Hamiltonian becomes 

H' = '^{p-A'f + U' + {p-A)V' + W'{p-A') 

= y(f+^j 



2 " ■ 2uq 



"^V + ^^T^^'- (3.12) 
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As a result, equation fl3.8p takes the form of equation for the harmonic oscillator: 

and can be solved, once again, by the traditional method of separation of the variables or by the 
factorization technique. 

4. Separation of Variables for a Shifted Harmonic Oscillator 
We shall refer to the case f l3.8p as one of a shifted linear harmonic oscillator. The Ansatz 

in the time-dependent Schrodinger equation results in the stationary Schrodinger equation 

H(f = E(f 

with the Hamiltonian (13.71) . The last equation, namely, 

- V? + X ip -\ [Zxip + ip) = — ip, 

with the help of the substitution 

fiXx^\ 
if = exp I u [x) 



2uj 



is reduced to the following equation 



— H — 5-x u = — u. 

^0 tUo 



4.1 



4.2 



4.3 



4.4 



4.5 



4.6 



UJ 

gives us the stationary Schrodinger equation for the simple harmonic oscillator 

v" + [2e ~i^)v = Q 

with e = E/uj, whose eigenfunctions are given in terms of the Hermite polynomials as follows 

and the corresponding eigenvalues are 

£n = n+^, En = uj(^+^ (n = 0, 1,2, ...) . 

Thus the normalized wave functions of our shifted oscillator (13. 8p are given by 

^„(x,t) = e--("+i/2)V„(a:), (4.10 

where 

^„ (x) = a exp ('-^) e-^'I'H^ (0 , i = x.[^ (4.11 
V J V ^0 



The change of the variable 




4.7 



4.8 



4.9 



and 



2 I 1 



\C\'= (4.12 
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in view of the orthogonality relation 

ipl{x)ip^{x) dx = Snm- (4.13) 

We use the star for complex conjugate. 

Solution of the initial value problem (12.151) can be found by the superposition principle in the 
form 

oo 

^lj{x,t) =J2cn^n{x,t), (4.14) 

n=0 

where 

oo 

^ (a;, 0) = X (x) = ^ c„ (x) (4.15) 

n=0 

and 



^*niy)xiy) dy (4.16) 

in view of the orthogonality property ( 14.130 . Substituting (14.160 into (14.14p and changing the order 
of the summation and integration, one gets 



^P{x,t) = / G{x,y,t)x{y) dy, (4.17) 

J — oo 

where the Green function is given as the eigenfunction expansion: 

oo 

G {x, y,t) = J2 e-^-("+^/2)V„ (x) < (y) ■ (4.18) 

n=0 

This infinite series is summable with the help of the Poisson kernel for the Hermite polynomials 
(Mehler's formula) [30]: 

i:^^^'--7Tbexp( ^-^-/f;^>\ M<i. (4.1, 

n=0 ■ V J- ' \ / 

The result is given, of course, by equation (13. 9p . 

5. The Factorization Method for Shifted Harmonic Oscillator 

It is worth applying the well-known factorization technique (see, for example, [2], [3], [1], [10] and 
[26] ) to the Hamiltonian (13. 7p . The corresponding ladder operators can be found in the forms 

d 

a = {a + iP) X + •y—, (5.1) 

d 

= (a — iP)x — 'j—, (5.2) 
ox 

where a, (3 and 7 are real numbers to be determined as follows. One gets 

aaV = (a' + /?')x2^ + (a-^/?)7^-2^/37^^-7'?t, (5.3) 

ox OX"^ 

a)ai) = (a'' + P^)x^ij-{a + iP)-f^-2ip-fx^--f^^, (5.4) 
^ ^ Ox Ox^ 
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whence 

(aa^ — a^a) ip = 2a'^il) (5-5) 

and 

1 (aat + ata) ^ = -7'^ + + - i^l (^x^ + ■ (5.6) 

The canonical commutation relation occurs and the Hamiltonian (13.71) takes the standard form: 

H {aa^ + a'^a) , (5.7) 

if 

2«7 = 1, u {a^ + P^) = u-f^ = ^uo, uPj = -^A. (5.8) 

The relation cUq = cu^ + A^, which defines the new oscillator frequency, holds. As a result, the explicit 
form of the annihilation and creation operators is given by 

v^a = [A J^ — , (5.9) 

sTla^ = (/Z+^),_ (5.10) 

The special case A = and uj = ujq gives a traditional form of these operators. 

The oscillator spectrum ( 14. 9p and the corresponding stationary wave functions ( 14. lip can be 
obtain now in a standard way by using the Heisenberg-Weyl algebra of the rasing and lowering 
operators. In addition, the n-dimensional oscillator wave functions form a basis of the irreducible 
unitary representation of the Lie algebra of the noncompact group SU (1, 1) corresponding to the 
discrete positive series P^; see [25], [28] and [33]. Our operators fl5.9l) - fl5.10l) allow us to extend 

these group-theoretical properties for the case of the shifted oscillators. We leave the details to the 
reader. 

6. Dynamics of Energy Related Expectation Values 
The expectation value of an operator A in quantum mechanics is given by the formula 

POO 

{A)= / ^p*ix,t) A{t)ip{x,t) dx, (6.1) 



where the wave function satisfies the time- dependent Schrodinger equation 



= Hij. (6.2) 
The time derivative of this expectation value can be written as 

^^^iA)^,{f)H^H-H'^)• (6.3) 

where is the Hermitian adjoint of the Hamiltonian operator H. Our formula is a simple extension 
of the well-known expression [21], [26], [32] to the case of a nonself- adjoint Hamiltonian. 
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We apply formula (16.31) to the Hamiltonian 

H = (p^ + x^) - Xpx, p = -i-^ (6.4) 
2 ^ ^ ox 

in equation (12. ip . A few examples will follow. In the case of the identity operator A = 1, one gets 

AH - H^A = X{xp- px) = iX (6.5) 
by the Heisenberg commutation relation 

[x,p] = xp — px = i. (6.6) 

As a result, 

^ "'' = A||V^f , (6.7) 



dt 

and time-evolution of the squared norm of the wave function for our model of the damped quantum 
oscillator is given by equation (12. 9p . 

In a similar fashion, ii A = H, then 

H"^ - H^H = {H-H^)H = iXH, (6.8) 

and 

^^{H) = X{H), {H) = {H),e'K (6.9) 

Moreover, 

j^{H-) = X{H-), (i/") = (i/")oe^* (n = 0,1,2,...), (6.10) 

which unifies the both of the previous cases. 

Now we choose A = p"^, A = x"^ and A = px + xp, respectively, in order to obtain the following 
system: 

A <^p2^ = SX {p^) - cuo (px + xp) , 

^{x^) = -X{x^) + ujq{px + xp) , (6.11) 



-f- (px + xp) = 2ujq ((p"^) — (x"^)) + X (px + xp) . 
dt \ \ / \ / / 



Indeed, 



^0 r 2 21 
[P 

?>iXp^ — iujQ {px + xp) 



p^H-H^p' = '^[p\x']+X[x,p^] (6.12) 



x^H — H^x^ = — [x'^,p'^~\ — Xx[x,p]x (6.13) 



2 

iujQ {px + xp) — iXx'^, 



and 



{px + xp) H — H'^ {px + xp) (6.14) 
2 



'^{[p,x'] + [x,p^]) 
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+"y {P [^^ p]p - X [x, p] x) 

+\ {{xpf - {px)^) 
2iuQ (p^ — x^) + i\ {px + xp) 



which results in (16.111) . 



The system can be solved explicitly, thus providing the complete dynamics of these expectation 
values. The eigenvalues are given by rg = A, r± = A±2iu; and the corresponding linearly independent 
eigenvectors are 









a^o = 1 








\ 2A j 





Ur 



(6.15) 



2ujq (A ± iuj) 



with the determinant 



ujq (A + iuj)^ (A — iu)' 



LUr 



(6.16) 



2A 2ijJq (A + iuj') Iujq (A — iuj) 
The general solution of the system (16.111) can be obtain in a complex form as follows 



ip') 
{px + xp) 

+C 6^^+^*^^^* 



Coe 



xt 



2A 



(6.17) 



LUr 



_|_ (J ^^{\-2iu))t 



(A — iuY 



UJn 



2u!o (A + iu) 



2uq (A — iu) 



where Cq and C± are constants. The corresponding solution of the initial value problem is given by 



ip') 
{px + xp) 



2c^2 



{uo ((p'>o + (x^) ) -X{px + xp)^) e^* ^0 

2A 



(6.18) 



1 / A 



2uj^ \ujq 



{px + xp)q + 



xe 



xt 



2uJoUJ 



xe 



xt 



(A^ — cu^) cos 2ujt — 2Xuj sin 2ujt 
cJq cos 2ujt 
2Xujq cos 2ujt — 2ujquj sin 2ujt 

{px + xp)r. - — (x^). 

2Aco' cos 2ujt + (A^ — ciJ^) sin 2ujt 
cUq sin 2ujt 
2ujquj cos 2ujt + 2Ati;o sin 2ujt 
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The mechanical energy operator E can be conveniently introduced as the Hamiltonian of our 
shifted linear harmonic oscillator fl3.7p : 

E = iJo = Y (p' + - ^ (px + xp) , (6.19) 

so that 

H = Ho + t^. (6.20) 



Then 



> - t(I«4(^=>) (6.21) 



X d , . 
-^j^ipx + xp) 

with the help of our system (16. lip . Therefore, 

^(E)=A(E), {E) = {E),e'' (6.22) 
for the expectation value of the mechanical energy of the damped oscillator under consideration. 

The case of the second Hamiltonian: 

^ = Y (p' + x') - Xxp = Ho- 2^, (6.23) 

which is the Hermitian adjoint of the Hamiltonian ( 16. 4p . is similar. Here 

jjn+l _ ^\^n =i^H- H^) H'' = X\p, x] H"" = -iXH"" 



and 



Moreover, 



±{H^) = -X{H^), {H^) = {H^),e-'' (n = 0, 1, 2, ...) . (6.24) 

p'^H - H^p^ = y [p^ x^] + Xp [x, p] p (6.25) 
= iXp'^ — iujQ {px + xp) , 



x'H-H^x' = ^[x\p^]+X[p,x'] (6.26) 



„2 rr tjt„2 _ ^ [^2 „21 

2 

= —3iXx'^ + iujQ {px + xp) , 

{px + xp) H — {px + xp) (6.27) 
= '^{[p,x'] + [x,p^]) 

+ ^ {p\x,p\p - X [x, p] x) 

—X {{xp)^ — (px)^) 
= 2iujQ (p^ — x^) — iX {px + xp) , 
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and the corresponding system has the form 

^ (/) = {P^) - ^0 {px + xp) , 

A^x^^ = -3A(x2)+c<;o(px + xp), (6.28) 
— {px + xp) = 2iUo {{p"^) — {x"^)) — A {px + xp) . 

UjL 

The change p ^ x, X ^ —A, ujq ~ujq transforms formally this system back into (16.111) . This 
observation allows us to obtain solution of the initial value problem from the previous solution given 
by fl6.18p . For the mechanical energy operator E introduced by equation f l6.19p one gets 

d_ 
dt 

with the help of (Km . 



-{E) = ~\{E), {E) = {E),e-'' (6.29) 



The case of a general variable quadratic Hamiltonian of the form 

H = a{t)p^ + h it) x^ + c it) px + d (t) xp, (6.30) 

where a{t) , b (t) , c{t) , d (t) are real- valued functions of time only, is considered in a similar fashion. 
One gets 

jjn+i _ ^\^n =(^H- H^) H'' = {c- d) [p, x\H'' = i{d- c) H"" (6.31) 

and 

m = l^) + {d{t)-c{t)){H-). (6.32) 



dt \ dt 

The cases n = and n = 1 result in 



and 

respectively. 
Moreover, 



(1) = (l)o exp (y^* {d (r) - c (r)) dr^ (6.33) 
JtW = (^^) + {d{t)-c{t)){H), (6.34) 



X 



p^H — H'^p^ = b [p"^, x^] + c [p^, x] + dp [p, x] p (6.35) 
= —i (3c + d)p^ — 2ib (px + xp) , 

H — H^x^ = a + ex a; + d (6.36) 
= i (3(i + c) + 2ia [px + xp) , 

{px + xp) H — [px + xp) (6.37) 
= a ([x,p^] +p[x,p\p) 
+b ( [p, x'^] + X [p, x] x) 
+ {c- d) ((px)^ - (xp)^) 
= Aiap^ — Aibx'^ — i{c — d) {px + xp) , 
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and the corresponding system has the form 

d 



{v^) = - (3c + d) {p^) - 2b {px + xp) , 

^ (x^) = (c + 3d) {x^) + 2a {px + xp) , (6.38) 

— (px + xp) = 4a (p^) — 46 (x^) + (d — c) (px + zp) . 
dt \ / \ / 

We have used the famihar identities 

[x,p\ = i, {xp)^ — ipx)"^ = i {px + xp) , (6.39) 
=2i{px + xp), [x,p^] =3ip'^, [x^,p] = Six'^ (6.40) 

once again. 

7. A Relation with the Classical Damped Oscillations 

Apphcation of formula (16.11) to the position x and momentum p operators allows to modify the 
Ehrenfest theorem [TT], [3^ for the models of damped oscillators under consideration. For the 
Hamiltonian (16. 4p one gets 

2 

pn — np = 

and 



xH — H^x = -^[x,p'^'\ = iuop, (7.1) 

pH — H^p = — [p, x^] + A [a;,p^] = —iuox + 2iXp (7.2) 
2 



^ {x) = ujo {p) , J^iP) = -^0 {x) + 2X {p) . (7.3) 



Elimination of the expectation value {p) from this system results in 



^ {x) - 2a| (x) + ul {x) = 0, (7.4) 



which is a classical equation of motion for a damped oscillator [3], [2U] . 
For the second Hamiltonian fl6.23p we obtain 



which gives 



in a similar fashion. 



^ (x) = ^0 {p) - 2A (x) , ^^P^ ^ ' ^^-^"^ 

^{x) + 2Xj^{x)+ujI{x) = (7.6) 



Finally, our model of the shifted harmonic oscillator (13. 8p . when the Hamiltonian is given by 
(I6.19p . results in 

^ (x) + {u;l - X') (x) = 0. (7.7) 

We leave the details to the reader. 
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8. The Third Model 



For the time- dependent Schrodinger equation with variable quadratic Hamiltonian: 



+ (8.1) 

where a = (a;o/2) e~^^*, b = {uJo/2)e'^^^ and c = d = 0, the characteristic equation takes the form 
fl2.2p with the same solution fl2.3p . The corresponding propagator has the form (11.21) with 

u cos ut — X sin ut n,., . . 

a(t) = e^^\ 8.2 

Pit) = ^e'\ (8.3) 

ujq sm ujt 

, , u cos ut + \ sin ut 

7W = ^ ■ 8.4 

Zujq sm oot 

This can be derived directly from equations (]1.2p -( fL8l) with the help of identity (12. 8p . We leave the 
details to the reader. It is worth noting that equation (18.11) can be obtain by introducing a variable 
unit of length x xe^^ in the Hamiltonian of the linear oscillator. 



9. Momentum Representation 



The time-dependent Schrodinger equations for the damped oscillators are also solved in the 
momentum representation. One can easily verify that under the Fourier transform our first Hamil- 
tonian (16.41) takes the form of the second Hamiltonian (16.231) with A —X and visa versa (see, 
for example, Ref. [9] for more details). Moreover, the inverses of the corresponding time evolution 
operators are obtained by the time reversal. Therefore, all identities of the commutative evolution 
diagram introduced in Ref. p] for the modified oscillators are also valid for the quantum damped 
oscillators under consideration. We leave further details to the reader. 
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